Spectral density in a nematic state of iron pnictides 
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Using cluster-perturbation theory, we calculate the spectral density A(k, cj) for a nematic phase 
of models describing pnictide superconductors, where very short-range magnetic correlations choose 
the ordering vector (7r,0) over the equivalent (0,7r) and thus, break the fourfold rotation symmetry 
of the underlying lattice without inducing long-range magnetic order. In excellent agreement with 
angle-resolved photo-emission spectroscopy (ARPES), we find that the yz bands at X move to 
higher energies. When onsite Coulomb repulsion brings the system close to a spin-density- wave 
(SDW) and renormalizes the band width by a factor of 2, even small anisotropic couplings of 
10 to 15 meV strongly distort the bands, splitting the formerly degenerate states at X and Y by 
^ 70 meV and shifting the yz states at X above the chemical potential. This similarity to the SDW 
bands is in excellent agreement with ARPES. An important difference to the SDW bands is that 
the yz bands still cross the Fermi level, again in agreement with experiment. We find that orbital 
weights near the Fermi surface provide a better characterization than overall orbital densities and 
orbital polarization. 

PACS numbers: 74.70.Xa,74.25. Jb,71.27.+a,71.10.Fd 



I. INTRODUCTION 

In recent years, iron-based superconductors have been 
intensely studied,^ ^ because of their high superconduct- 
ing transition temperatures. As in the cuprates, antifer- 
romagnetic (AFM) order is present in the phase diagram 
and phonons are not believed to be strong enough to ex- 
plain the high transition temperatures.^ However, in the 
pnictides the AFM phase is a metallic spin-density wave 
(SDW) rather than a system of localized Heisenberg spins 
as is the case in the cuprates. At temperatures slightly 
above the transition to the SDW with ordering vector 
(tt, 0) or (0,7r), many weakly doped compounds show an 
orthorhombic phase without long-range magnetic order, 
but with broken rotational symmetry. This phase has 
slightly different lattice constants along the in-plane iron- 
iron bondsj^ but the anisotropy that develops in elec- 
tronic observables such as resistivit}'^ ^ or angle-resolved 
photoemission spectroscopy (ARPES)P'^^ of detwinned 
samples appears considerably more pronounced. 

A number of competing scenarios have been proposed 
for this phase and can be broadly categorized as "mag- 
netism", "orbital", and "lattice" driven. In the first case, 
the symmetry between equivalent magnetic ordering vec- 
tors ('7r,0) and (0,7r) is broken and the system chooses 
one of them wi thout immediately establishing long-range 
magnetic order .HH^ In the second picture, it is the degen- 
eracy between two d orbitals of the iron ion, the xz and yz 
states providing the greatest contribution to the states at 
the Fermi surface (FS), that is spontaneously broken; 
the resulting orbital occupation then determines the ef- 
fective magnetic exchange constants that generate the 
SDW order. Both pictures were first discussed in insu- 
lating spin and spin-orbital models and have since been 
generalized to take into account electron itineracy. Stud- 



ies in several models have shown that nematic phases can 
indeed develop between structural and magnetic transi- 
tion temperatures.^^ 

While a definite answer about the driving mechanism 
ma y b e hard to nail down, as spin, orbital, and lat- 
grees of freedom are naturally coupled and in- 
teract with each other, one may nevertheless try to iden- 
tify the dominant ingredient (s). To this end, it is instruc- 
tive to establish how each type of symmetry breaking 
manifests itself in observables. If the symmetry breaking 
is assumed to mostly concern the xz and yz orbitals, one 
can introduce it explicitly by adding a phenomenologi- 
cal energy splitting between the orbitals and evaluating 
its impact on observables such as the optical conductiv- 
ity or the spectral density. These signatures were found 
to qualitatively agree with experiments,^^ where states 
of yz character are found to be higher in energy than 
those of xz character in several different pnictide com- 
pounds from the two structurally slightly different "111" 
and "122" families. On the other hand, ARPES data 
taken above the Neel temperature have alternatively been 
interpreted in terms of "band folding" due to magnetic 
ordei^ or emphasising the coupling between magnetic or- 
der and the orbital states near the Fermi level. Short- 
range magnetic ordei^^ and the spin-nematic scenaricP^ 
likewise reproduce the anisotropic conductivity, and the 
latter has been argued to lead to an effective orbital split- 
ting.^^ However, a direct calculation of the spectral den- 
sity in a nematic phase is so far lacking. 

We use here a method that combines real and mo- 
mentum space, cluster-perturbation theory (CPT),I^2^ 
to calculate the one-particle spectral density 74(k,a;) 
for a spin-nematic phase where rotational symmetry 
is broken via (very) short-range spin correlations that 
are AFM in x and ferromagnetic (FM) in y direction 
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[corresponding to ordering vector (7r,0)], but without 
long-range magnetic correlations beyond second neigh- 
bors. The obtained spectral density reproduces the 
momentum-dependence of the band shifting observed 
in ARPES.^^^^ If onsite interactions bring the system 
close to the SDW transition, a small phenomenologi- 
cal magnetic anisotropy leads to large anisotropics in 
A(k, cj). Thus, we can theoretically describe the aston- 
ishing ARPES result that the overall band shifts charac- 
terizing 74(k, u) in the SDW phase are nearly fully devel- 
oped already above the Neel temper at ure.l^ 

In order to be able to solve the problem exactly on 
a small cluster, we use variants of models with thre^^ 
and fouil^ orbitals; models and method are introduced 
and discussed in Sec. HIl Section IIII Al discusses the 
anisotropic band shifts induced by (strong) anisotropic 
magnetic couplings in the non-interacting models; in 
Sec. |IIIB[ we show that in the presence of onsite inter- 
actions and near the SDW transition, smaller magnetic 
anisotropics have a large impact. In Sec. 
are summarized and discussed. 
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II. METHOD AND MODEL 

The aim of the paper is to calculate the spectral density 
yl(k, cj) in a phase where (short-range) magnetic corre- 
lations break the fourfold symmetry of the lattice, but 
without long-range magnetic order. The latter require- 
ment prevents us from carrying out our calculations di- 
rectly in momentum space, as it has been done for the 
paramagnetic and AFM phases. As an alternativ e ap- 
proach, we choose here cluster perturbation theory.^^^ 
In this method, the ground state and one-particle Green's 
function are evaluated (almost) exactly (with Lanczos 
exact diagonalization) for a fully interacting quantum 
model on a small cluster, and hoppings between clusters 
are treated in perturbation theory [for an illustration see 
l(b)| . Apart from the limit of small inter- 



Figs. 1(a) and 



cluster hoppings, this approximation also becomes exact 
in the opposite limit of vanishing interactions, as can be 
seen by considering that it amounts to replacing the self 
energy of the full system by that of the small cluster.!^ 
Long-range order can be treat ed wit h the related varia- 
tional cluster approach (VGA) pHEll as it has been done 
for a two-orbital model for pnictides.'^^^ 

The biggest drawback of the VGA is that correlations 
are only included exactly within the small cluster, while 
longer-range effects are treated at a mean-field level. For 
nematic phases with at most short-range order, this lim- 
itation turns into a huge advantage: We can break the 
symmetry between the x and y directions locally on the 
small cluster, see below, but without imposing long-range 
order by a symmetry-breaking field. If the small cluster 
is, e.g., an AFM coupled dimer, its groundstate is thus 
still given by a singlet, i.e., a superposition of "up-down" 
and "down-up", which removes long-range correlations. 

When using a dimer as the directly solved cluster, we 



FIG. 1. (Color online) Schematic illustration of the cluster de- 
composition into (a) four-site and (b) three-site clusters used 
for the three- and four-orbital models. Ground-state ener- 
gies and Green's functions of the small cluster - as connected 
by thick solid lines - are obtained by exact diagonalization. 
Clusters are then connected in GPT along the thinner dashed 
bonds. Within the cluster, AFM (FM) Heisenberg exchange 
acts between electrons along the x- (y-) bond, (c) Spec- 
tral density A(k, cj) of the non-interacting three-orbital model 
Eq. Q . Solid lines indicate the bands in terms of the pseudo- 
crystal momentum k, shading the spectral weight in terms 
of "real" momentum k, the difference is that weight with xy 
character shifts by (tv^tv). [Note that along (0, 0)-(7r, tt), the 
bands with dominant xz/yz character contain these two or- 
bitals with identical weight, even though the yz character - 
here drawn on top - dominates the figures.] (d) A(k, cj) of 
the non-interacting four-orbital model. Dashed lines indicate 
the results for the model that is obtained by removin^^ the 
— orbital from the five-orbital model of Ref. 1251 Solid 
lines are for the model used here, where the x^ — y^ orbital is 
then somewhat removed from the Fermi level by changing t^^ 
from -0.02 to tH = 0.03 and 63 from -0.22 to -0.12 (nota- 
tion as in Ref. 1^5)) . As in (c), shading indicates the spectral 
weight for the "real" momentum instead of the pseudo-crystal 
momentum. In the online version, red refers to xz, blue to 
yx, and green to all other orbitals. In all spectra, peaks are 



broadened by a Lorentzian 6/{{(jJ 



' + S^) with S = 0.05 



except for Fig. 1(d), where 6 = 0.025. All energies are in eV. 



find instabilities, i.e., poles of the one-particle Green's 
function that are on the wrong side of the chemical po- 
tential. While this does not necessarily invalidate the 
results (which are in fact similar to the more stable re- 
sults described below), it may indicate that the self en- 
ergy of a dimer differs too strongly from that of a large 
two-dimensional system to provide a reliable approxima- 
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tion. In order to be able to use three- (four-) site clus- 
ters instead, which lead to stable results, we restrict the 
Hamiltonian to the four (three) orbitals that contribute 
most of the weight at the FS. The results presented here 
were obtained with the cluster decompositions shown in 
Figs. |l(a)| and |l(b)[ but equivalent results were found 
for the three-orbital model when using a "brick-wall" ar- 
rangement of 2 X 2 clusters instead of the "columns" in 

The momentum-dependent tight-binding Hamiltonian 
in orbital space can be written as 

where d- (dl ) annihilates (creates) an electron 

with pseudo-crystal momentum k and spin a in orbital 
ly. The three-orbital model used here is based on the 
model of Ref. 23, but a few longer-range hoppings were 
added to provide a better fit of the bands near the FS, 
because the original three-orbital model has magnetic in- 
stabilities too far from (tt, 0)/(0, tt).!^ The T^'^(k) give 
the hoppings between orbitals /i and v and are 

2^11/22 _ 2^2/1 COS kx + 2^1/2 COS ky + 4^3 COS kx COS ky 

± 2^11 (cos 2kx — cos 2ky) + 4ti2 cos 2kx cos 2ky^ 

(2) 

T^^ = Axy + 2^5 (cos kx + cos ky) + Ate cos kx cos ky 
+ 2^9 (cos 2kx + cos 2ky) 

+ 4tio(cos 2kx cos ky + cos kx cos 2ky)^ (3) 

2^13 ^ fsi ^ 2itY sin kx + 4it8 sin kx cos A:^^, (5) 
2^23 ^ j-32 ^ 2it7 sin ky + 4it8 sin fc^ cos kx, (6) 

where a bar denotes the complex conjugate. Hopping 
parameters are h = -0.08, t2 = 0.1825, ts = 0.08375, 
t4 = -0.03, h = 0.15, te = 0.15, = -0.12, = 
-t7/2 = 0.06, tio = -0.024, til = -0.01, ti2 = 0.0275, 
Aa,^; = 0.75, /i = 0.4745; Fig. Jl] shows the uncorre- 
lated tight-binding bands. The four-orbital model was 
obtained by removing the Sz"^ — orbital^^ from the 
five-orbital model of Ref. 25 and slightly changing onsite 
energy and third-neighbor hopping of the — y'^ orbital 
to alleviate the fact that removing the 3z^ — orbital 



moves it too close to the Fermi level, see Fig. |l(d) 



In 

principle, hoppings can be extended to three dimensions 
and parameters could be fitted to model specific com- 
pounds, at least in the more detailed four-orbital model. 
The features we aim to study here - an anisotropy be- 
tween the X and Y points - have been experimentally 
observed in different compounds, and we are going to see 
that both the three- and the four-orbital models lead to 
similar results despite their somewhat different disper- 
sions, suggesting that fine-tuning of the kinetic energy is 
not crucial. 

We use a unit cell with one iron atom to distinguish be- 
tween momenta (tt, 0) and (0, tt), which would both map 



to (tt, tt) for a two-iron unit c ell. D ue to an internal sym- 
metry of the two- iron unit cell j^^l^^4 t is always possible to 
use a one- iron unit cell for tight-binding models restricted 
to an Fe-As plane. However, the xz and yz orbitals with 
momentum k couple to the other orbitals at momentum 
k+ (tt, tt). Thus, one writes the tight-binding Hamiltoni- 
ans in terms of a pseudo-crystal momentum k, which is 
k = k for xz/yz and k = k+(7r, tt) for xy / x'^ —y'^ / 3z'^ —r'^ . 
In real space, such a notation corresponds to a local gauge 
transformation, where replacing, e.g., the xy[ orbital at 
site i = {ix^iy) by (— l)^*^+*^^^xyi (and analogously for 
x'^ — y'^ and 3z'^ — r^) leads to a translationally invariant 
Hamiltonian with a one-iron unit cell. For comparison 
with ARPES experiments, however, this gauge transfor- 
mation has to be undone, which implies that spectral 
weight at k with orbital character xy, x^ — y^ or Zz^ — 
is plotted at k = k + (tt, 7r).C^ 

In order to study a nematic phase, the four- fold lat- 
tice symmetry is explicitly broken by introducing a 
phenomenological Heisenberg interaction, which couples 
electrons in all orbitals, 



Heis 



(7) 



(iJ)ll^/y 



where /i, u denote orbitals and nearest-neighbor 
(NN) bonds. For Janis > 0, the coupling is AFM 
(FM) along the x {y) direction. The electron-spin op- 
erators are given by Si^ = ^Y.ss' ^Ls^ss'^ius'^ where 
a = {a^ ^ (j^) is the vector of Pauli matrices. These in- 
teractions act only within the small cluster that is solved 
exactly. We are here not going to investigate the origin of 
such a breaking of rotational sym metry, whi ch has been 
shown to occur in several mo dels ,^^^^'1^^ but we will 
study its impact on the system. We find that when the 
system is close to the spin-density wave, very small val- 
ues of Janis trigger highly anisotropic band distortions, 
suggesting that short-range correlations, as observed in 
a spin-fermion model, indeed favor such a symmetry 

breaking. 

When onsite interactions^^'^^" are taken into account, 
the same values of intra-orbital Coulomb repulsion /7, 
inter-orbital repulsion Hund's rule coupling J and 
pair hopping J' = J were used for all orbitals, along 
with the standard relation U' = U — 2 J, giving 



^int = U ^ni,a,t^i,a,i + {U' - J/2) ^ ni,c,ni,/3 

i,Q! i,Q;</3 
i,Q!</3 



(8) 



i,a</3 



where denote the orbital and S\^cx (^i,a) is the spin 
(electronic density) in orbital a at site i. While the pa- 
rameters relating to the xy and x'^ — y^ orbitals can in 
principle be slightly different from each other and the 
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diflFerences between X and Y are large, changes around 
r = (0,0) are far less pronounced. Th e corr esponding 
orbital-resolved FS can be seen in Fig. 
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FIG. 2. (Color online) (a) Spectral density A(k,(jj) and 
(b) Fermi surface of the three-orbital model (four-site clus- 
ter) with parameters as given for Fig. l(c)| and an explicit 
symmetry-breaking Janis = 0.5 eV, see Eq. but with- 
out Coulomb repulsion and Hund's rule coupling. Shadings 
are for "real" momentum, lines indicate the non-interacting 
model in pseudo-crystal momentum k. 



xz/yz doublet, symmetric interactions were chosen for 
simplicity. 



III. RESULTS 



A. Band anisotropy in the three- and four-orbital 
models 



In order to study the effects of phenomenological short 
range magnetic correlations, the Hamiltonian given by 
Eqs. ([T]) and ([7|, was initially treated with the VGA on a 
four-site cluster, with AFM interactions along x and FM 
ones along y but without onsite Coulomb and Hund inter- 
actions. A fictitious chemical potential was optimized as 
a variational parameter, but did not have a large impact 
on the results. No tendencies towards long-range order 
were found, which agrees with expectations: Since the 
AFM Heisenberg interaction only acts within the cluster, 
it favors a total cluster spin of 5'tot = 0. In the large sys- 
tem, consisting of many noninteracting clusters, there is 
no magnetic order. Rather large Janis ^ 0.3 eV has to be 
chosen to induce appreciable signatures of the anisotropy, 
which is a very large energy scale compared to the other 
parameters of the Hamiltonian. The reason is that the 
non-interacting model with four electrons per site does 
not contain any net unpaired spins that can directly be 
coupled by a Heisenberg interaction; the interaction first 
needs to be strong enough to induce a local spin. 

The spectral density for Janis = 0.5 eV is shown in 
Fig. |2(a)[ Apart from the fact that the Heisenberg inter- 
actions make the spectrum more incoherent, the bands 
are most strongly modified near X = (7r,0), which cor- 
responds to the ordering vector that would be favored 
by the NN AFM interaction along x. One clearly sees 
that the yz states around X are moved to higher en- 
ergies, while the xz states at F = (0, tt) are shifted to 
slightly lower energies in agreement with experiments. 
The energy shifts are momentum dependent: While the 



Like the 

spectral density, it shows some features that are similar to 
those resulting from band folding in a (tt, 0) SDW; for ex- 
ample, the xz electron pocket at Y has a "mirror pocket" 
at M = (tt, tt) = 7" + (tt, 0). However, the FS is still qual- 
itatively different from the FS found in the long-range 
ordered SDW, where folding leads to additional features 
and largely suppresses the yz weight, which dominates 
the hole pockets here. Such differences related to long- 
range order are consistent with ARPE S expe riments in 
NaFeAs,^ see also the discussion in Sec. IIIB 



The same behavior as in the three-orbital model is seen 
for the four-orbital case, see Fig. [sj where A(k, cj) is 
shown for increasing Janis = 0.2, 0.3, 0.4 eV. In the 
last case, the splitting between the states at X and Y is 
^150 meV. Taking into account that the overall band 
width has to be renormalized by a factor of 2-3, this is 
consistent with the order of magnitude of the 60 meV 
splitting reported for Ba(Fei_ajCoaj)2As2f^ This can be 
compared to an explicit orbital splitting, similar to the 
mechanism proposed in Ref. [181 The splitting can be 
written as A = {nyz—rixz)/'^^ where Uxz i^yz) is the den- 
sity in the xz (yz) orbital, and was set to A = 0.15 eV, 
which approximately reproduces the energy difference be- 
tween t he X and Y points indicated by the dashed lines 
in Fig. 3(c) A momentum-independent splitting large 



enough to reproduce the energy differences between the 
X and Y points substantially distorts the features near 
the F point as well. While the unoccupied states above 
the chemical potential at F are not easily accessible in 
ARPES, available data on the bands defining the hole 
pockets appear more consistent with the slighter changes 
caused by momentum-dependent shifts of the nematic 
scenario, especially for cases with a large splitting be- 
tween the features near X/Y", where the bands near F 
would be very strongly distorted by a rigid shift. ^ "^^ 

Total orbital densities do not turn out to be a reliable 
way to characterize the impact of the nematic order on 
states near the Fermi energy. Densities in the xz and yz 
orbitals differ only slightly in the four-orbital model with 
^xz — ^yz ^ 0.02 for Janis = 0.4 cV. This value is not 
strongly affected by 4% hole or electron doping, in con- 
trast to a proposed sign change for hole doping^^ and it 
is broadly consistent with the small orbital polarizations 
found in mean-field analyses for the SDW state. In 
the three-orbital model, the orbital polarization is even 
opposite with Uxz — ^yz ^ —0.1, because spectral weight 
with yz character is transferred bel ow th e Fermi level [see 
the density of states shown in Fig. 4(a)|, and in contrast 
to the four-orbital model [see Fig. 4(b) , this weight is 
not balanced by xz states further away from ji. Never- 
theless, the band reconstruction near the Fermi level and 
the band anisotropy are very similar in the two mod- 
els. AFM correlations along x always bring the yz states 
around X closer to the Fermi level, even when the to- 
tal orbital densities satisfy Uyz > Uxz^ in contrast to 
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FIG. 3. (Color online) Spectral density A(k, cj) of the four- 
orbital model (three-site cluster), see Fig. 1(d) and an 
increasing explicit symmetry-breaking term Eq. (T]) of (a) 
Janis = 0.2 eV, (b) Janis = 0-3 cV, and (c) Janis = 0-4 cV. 
Coulomb repulsion and Hund's rule coupling are not included. 
Shadings are for "real" momentum, solid lines indicate the 
non-interacting model in pseudo-crystal momentum k. In (c), 
dashed lines are for a non-interacting model with an energy 
difference A = 0.15 eV between the xz and yz orbit als, which 
was fitted to approximately reproduce the difference between 
the X and Y points. 



a naive expectation that the yz bands should be low- 
ered in energy in this case. When onsite interactions 
bring the thr ee-orb ital model closer to the SDW transi- 
tion (see Sec. Ill B| below), the orbital densities become 
a,,, ^ -0.012 for /7 = 1 eV. To- 



almost equal with — riyz 
tal orbital densities can determine magnetic properties 
via the Goodenough-Kanamori rules in Mott insulators, 
which do not have a FS. On the other hand, since the 
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FIG. 4. (Color online) Density of states for (a) the three- 
orbital model with Janis = 0.5 eV and (b) the four-orbital 
model and Janis = 0.4 eV. U = Jnund = in both cases. 



pnictides are more metallic with strongly hybridized or- 
bit als, a more consistent and clearer picture can here be 
obtained if one concentrates on spectral weight near the 
Fermi level as will be discussed below. 



B. Impact of onsite Coulomb interaction 

In this subsection the impact of onsite interactions will 
be investigated. The full Eq. ([8| including spin-flip and 
pair-hopping terms can easily be included in the VGA. 
Interaction strengths were chosen below the critical val- 
ues for the onset of long-range order because we want to 
focus on short-range correlations here. As can be seen in 
Fig. |5] for the two models considered here, lower values of 
-^anis ^ 0.2 is now sufficient to induce substantial asym- 
metries, in contrast to the larger Janis ^ 0-4 to 0.5 eV 
needed for the noninteracting models. Onsite interac- 
tions favor local magnetic moments, even in the absence 
of long-range order, that can then be coupled even by 
weaker Janis- 

Finally, we study the three-orbital model very close 
to the SDW transition by setting U = 1.02 eV. In a 
mean-field treatment as used in Ref. [23l one finds an 
SDW with long-range magnetic order, but the optimal 
VGA solution does not yet show long-range order due 
to the presence of quantum fluctuations. However, the 
system is so close to a magnetically ordered state that 
very small Janis ^ 0.015 eV = 15 meV already intro- 
duces strong short-range order and corresponding band 
anisotropics. Several occupied low-energy bands [e.g. 
between F = (0,0) and Y = (0,7r) as well as around 
M — (tt, tt)] in the spectral density, which is shown in 
Fig. 5(c)[ have energies reduced by a factor of ~ 2, con- 
sistent with the renormalization factor ^2-3 needed to 
reconcile density- functional bands with ARPES. Bands 
above the Fermi level do not have reduced widths. This 
asymmetric impact of correlations is in agreement with 
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FIG. 5. (Color online) Spectral density with anisotropic 
short-range magnetic order and onsite interactions, (a) For 
the four-orbital model and U = 0.3 eV, Jnund = 0.075 eV 
and Janis = 0.2 eV. (b) For the three-orbital model with 
U = 0.6 eV (Jnund = 0.15 eV) and Janis = 0.2 eV and (c) 
[/ = 1.02 eV (Jnund = 0.255 eV) and Janis = 0.015 eV. For 
these last values of U and Jnund, the three-orbital model is 
very close to the SDW. (d) FS corresponding to the param- 
eters in (c), it captures spectral weight within 1 meV of the 
Fermi level; broadening of the spectral weight is consistent 
with (c). Onsite interactions are lower for the four-orbital 
model, because it is at half filling and Hund's rule thus, moves 
it closer to a Mott transition, while it partly compensates U 
away from half filling, as in the three-orbital model. Shadings 
are for "real" momentum; lines indicate the non-interacting 
model in pseudo-crystal momentum k. 



dynamical mean- field studies.*^ In addition to the renor- 
malization, Janis = 15 meV induces an energy splitting 
of ~ 70 meV between the X and Y points. In fact, the 
band at X has moved slightly above the chemical poten- 
tial, as expected for the SDW phase. The fact, that this 
happens even in the absence of long-range order, is in ex- 
cellent agreement with recent ARPES data for NaFeAs, 



where it was likewise found that the overall band posi- 
tions at X and Y nearly reach their "SDW values" above 
the Neel temperature.^ 

Nevertheless, the corresponding FS, see Fig. |5(d)[ 
clearly shows important differences to that of the SDW 
state: As the yz states cross the chemical potential here 
with a rather low Fermi velocity (leading to an elonga- 
tion of the hole pocket at F along the x-direction), they 
contribute substantial weight to the FS. In fact, both of 
the strong features along the T-X line are of yz char- 
acter. In the SDW phase, in contrast, the yz orbital 
dominates the AFM order parameter and is thus mostly 
gapped out Related effects have likewise been observed 
in ARPES, wher e these yz bands open gaps at the Neel 
temperature's 



IV. SUMMARY AND CONCLUSIONS 

The variational cluster approach was used to study 
the spectral density of a nematic phase in three- and 
four-orbital models for iron-based superconductors. We 
found that the method is well suited for problems involv- 
ing short-range correlations without long-range (mag- 
netic) order. The correlations considered in this study 
were extremely short-range, going only over NN sites, 
the minimum to break rotational invariance. While this 
is a somewhat extreme scenario, it has been argued that 
magnetic correlations that are effective only on a very 
short range lead to the linear temperature dependence of 
the magnetic susceptibility at high temperatures.!^ Nu- 
clear quadrupole resonance^^ indicates that there are As 
ions seeing different electronic surroundings in the "un- 
derdoped regime", which would be in agreement with 
the present scenario of As ions involved in "magnetic" 
vs. "non-magnetic" bonds. 

When symmetry between the x- and y directions is 
broken by a phenomenological magnetic interaction that 
is AFM in the x direction, the bands with yz charac- 
ter around momentum X = (tt, 0) move to higher en- 
ergies, i.e., closer to the Fermi level. This is in agree- 
ment with ARPES on detwinned samples above the mag- 
netic transition temperature, in both the "122" com- 
pound Ba(F ei_a;Coa;)2As2P and the "111" compound 
NaFeAs.'^Eni xhe latter is not expected to have surface 
state^^ that might complicate the analysis of ARPES 
in 122 compounds. The changes in the band structure 
due to the nematic order not only depend on the orbital, 
but also on momentum. Changes around the F = (0,0) 
point are far less pronounced than differences between 
X and r, again in agreement with ARPES. ISKni Total 
orbital densities and their difference are model depen- 
dent and not a reliable predictor of reconstructions of 
low-energy states. However, the orbital-resolved spectral 
weight and the bands near the Fermi level are affected in 
the same way both in a three- and a four-orbital model, 
with and without onsite interactions indicating that they 
are more universal and less dependent on details of the 
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model Hamiltonian. In agreement with previous find- 
ings on the orbital polarization of t he FS in the SDW 
phas^^ and on transport propertiesf^^E^ this suggests 
that total (orbital) densities are here less important than 
in Mott insulators, as the metallic character of the pnic- 
tides makes states near the Fermi level far more impor- 
tant than those further away. 

When onsite interactions are strong enough to bring 
the system close to the SDW transition, very small 
anisotropic couplings can deform the bands until their 
broad features resemble bands in the SDW regime, i.e., 
bands are renormalized by a factor of ~ 2 and the yz 
states at X move above the chemical potential, as seen 
in ARPES on NaFeAs just above the Neel temperature.^ 



Nevertheless, the Fermi surface still differs from that of 
a state with full long-range magnetic order, where the 
yz states are mostly gapped out j^^ * ^^ * again in agreement 
with ARPES.I^ 
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